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PREFACE 

The  impetus  for  an  investigation  of  the  world-wide  incidence  of  daily 
freezing  conditions,  especially  diurnal  freeze-thaw  cycles,  was  generated  by 
requests  from  several  Army  agencies,  including  the  Test  and  Evaluation  Command 
(TECOM).  Freeze-thaw  cycles  ha*e  numerous  applications  of  military 
significance,  particularly  with  regard  to  off-road  mobility,  vulnerability  of 
roofing,  and  equipment  malfunction.  Climatological  studies  of  this  nature 
provide  weather  support  for  Field  Army  tactical  operations. 

This  work  was  accomplished  under  Project  4A1 6 1 102AT2U,  Task  C,  Work  Unit 
001,  "Relationship  between  Environmental  Factors  and  Materiel  Design 
Problems. " 

Appreciation  is  extended  to  the  Environmental  Technical  Applications 
Center,  Asheville,  N.C.,  for  providing  processed  data  for  many  stations  in  the 
Northern  Hemisphere;  to  Mark  Schroeder,  Mary  McClarnon,  Cedric  Key,  and  Lee 
Morkes,  U.S.  Army  Engineer  Topographic  Laboratories  (ETL),  for  computer  and 
drafting  assistance;  and  to  my  colleagues  in  the  Battlefield  Environmental 
Effects  Group,  ETL,  for  their  helpful  suggestions. 

The  work  was  performed  under  the  supervision  of  D.  W.  Dery,  Chief, 
Battlefield  Environmental  Effects  Group  (BEEG),  R.  J.  Orsinger,  Chief,  Land 
Combat  Systems  Division,  and  W.  E.  Boge,  Director,  Geographical  Sciences 
Laboratory. 


Col  Edward  K.  Wintz,  CE  was  Commander  and  Director  and  Mr.  Walter  E. 
Boge,  was  Technical  Director  of  the  U.S.  Army  Engineer  Topographic 
Laboratories  during  the  report  preparation. 
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DIURNAL  FREEZE-THAW  FREQUENCIES 
IN  SELECTED  REGIONS  OF  THE  HIGH  LATITUDES 


INTRODUCTION 

The  impact  of  diurnal  freeze-thaw  cycles  on  soils,  rocks,  roadbeds,  and 
construction  materials  is  of  much  concern  to  the  agriculturist,  hydrologist, 
geologist,  and  transportation  or  construction  engineer.  As  noted  by  Troll, 
such  alternation  of  freezing  and  thawing  affects  the  particular  size  structure 
of  the  soil,  causes  erosion,  runoff,  or  flooding,  and  reduces  the  bearing 
strength  of  the  surface  layer.1 2  The  vulnerability  of  the  ground  depends, 
among  other  variables,  on  the  amount  of  available  moisture  and  whether  the 
ground  is  bare  or  is  covered,  as  with  vegetation,  ice,  or  snow. 

The  surface  covering  affects  the  albedo,  which  in  turn  is  largely 
responsible  for  local  differences  in  extreme  temperatures,  i.e.,the  daily 
maximum  and  minimum,  the  determining  factors  of  daily  freeze-thaw.  Daily 
freezing  conditions  may  be  defined  as  consisting  of  frost  days  (min  <0°C), 
ice  davs  (max  <  0°C) ,  and  freeze-thaw  days  (min  <  0°C,  max  >0°C).  The 
interrelationship  among  these  three  variables  forms  the  basis  of  this  study. 
The  principal  question  is:  What  incidence  of  diural  freeze-thaw  cycles  may  be 
expected  at  a  given  site  per  given  interval  of  time? 

Much  of  the  information  in  this  report  was  presented  at  the  Fourth 
International  Conference  on  Permafrost,  held  18  through  22,  July  1983  at 
Fairbanks,  Alaska.  The  conference  paper  was  entitled  "Diurnal  Freeze-Thaw 
Frequencies  in  the  High  Latitudes:  A  Climatological  Guide." 

Information  is  provided  on  daily  freezing  conditions  in  the  permafrost 
and  contiguous  regions  of  Alaska,  Eastern  Siberia,  Iceland,  and  Greenland. 

Also  included  are  several  German  stations,  representing  a  range  of  elevations, 
in  order  to  compare  the  effects  of  altitude  with  those  of  latitude.  As  in  an 
earlier  study, ^  a  number  of  guides  are  offered  for  estimating  the  respective 
frequencies  of  frost  days,  ice  days,  and  freeze-thaw  days  per  month  or  year 
for  a  given  site  from  routine  climatological  parameters. 


lC.  Troll.  "Structure  Soils.  Solifiuction,  and  Frost  Climates  ol  the  Earth,"  S1PRE  Translation,  1958. 

2 

R.L  Wexlcr.  “A  General  Climatological  Guide  to  Daily  Freezing  Conditions:  Frost  Days,  Ice  Davs,  and  Freeze-Thaw  Days” 
U.S.  Army  Corps  of  Engineers,  Fort  Delvoir,  Virginia.  ETL-0287,  AD-A116  771,  1982.' 
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BACKGROUND 


The  geographical  distribution  of  the  annual  number  of  diurnal  freeze-thaw 
cycles  has  been  determined  for  various  countries  or  sections:  United  States, 
Canada,  Poland,  Japan,  the  Arctic,  Europe,  and  the  U.S.S.R.1  Annual  or 
monthly  frequencies  of  frost  days,  ice  days,  and  freeze-thaw  days  have  been 
correlated  respectively  with  mean  daily  minimum  temperatures,  mean  daily 
maximum  temperatures,  or  a  combination  of  both.  * 


DATA 

Freauencies  of  diurnal  freeze-thaw  cycles  are  not  readily  available.  Many 
climatic  sunmaries  list  frost  days,  but  few  list  ice  days.  In  the  past 
information  on  freeze-thaw  cycles  has  been  simulated  (Hastings,3  Visher’14). 
In  this  report,  all  the  analyses  are  based  on  observations  of  daily  maximum 
and  minimum  temperatures,  the  data  for  which  were  obtained  from  a  variety  of 
sources,  mginly  the  U.S.  Department  of  Commerce,^  the  Danske  Meteorologisk 
Institute,  and  the  U.S.  Air  Force  Environmental  Technical  Applications 
Center.'  All  the  temperatures  referred  to  in  this  report  were  from  standard 
weather  shelters  at  1.5  to  1.8  meters  above  ground  (screen  height). 


Refer  to  the  References  for  information  on  diurnal  freeze-thaw  cycles* 
Fraser,  1959;  Hastings,  1961;  Hershfield,  1972;  Pelko,  1970;  Russell,  1942; 
Rhit.ara,  1970;  Visher,  19^5 ;  Wevi«n,  1982;  end  Williams,  1964 

•  »  _ 

Refer  to  the  References  for  information  on  frequencies:  Fraser,  1959* 
Hershfield,  1972;  Shitara,  1970;  and  Wexler,  1982 


^A.D-  Hastings.  "Atlas  of  Arctic  Environment  "  Headquarters  Quartermaster  Research  and  Engineering  Command,  Natick, 
Massachusetts.  11)61. 

^S.S  Visher.  "Climatic  Maps  of  Geologic  Interest."  Bulletin  of  the  Gcnlugic  Society  of  America,  vol.  56,  pp.  713-736,  1915 

■’U.S.  Department  of  Commerce.  "Local  Climatological  Data  Annual  Summary  with  Comparative  Data  (stations  in  Alaska)." 
1980, 

^Danske  Metcorologiskc  Institute,  Meteorologisk  Arbog.Part  2,  Greenland,  1947-1965. 
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FROST  DAYS  AND  ICE  DAYS 


For  a  network  of  stations  in  a  given  region,  summarized  observations  per 
month  and  year  were  obtained  for  the  mean  daily  minimum  temperature,  N;  the 
mean  daily  maximum  temperature,  M;  the  number  of  frost  days,  F;  and  the  number 
of  ice  days,  I.  From  these  observations,  simple  linear  regression  models  were 
determined  such  that  F  may  be  derived  from  N,  and  I  may  be  derived  from  M  for 
any  site  within  the  specified  region  per  given  interval  of  time.  Figure  1 
gives  examples  of  these  regression  plots  for  annual  data  for  Alaska,  Eastern 
Siberia,  Iceland,  and  Greenland.  Figure  P  contains  similar  plots  for  monthly 
data  for  May  (Greenland)  and  April  (the  other  areas).  All  the  days  of  the 
month  were  usually  below  freezing  if  the  mean  maximum  temperature  was  <-6°C 
and  above  freezing  if  the  mean  minimum  temperature  was  >6°C.  The 
corresponding  regression  equations  are 


F  =  a1  -t-  bnN  (1) 

and 

I  =  &2  +  b£M  (2) 


The  parameters  a_  and  b_  depend  on  the  data.  They  serve  as  constants  for 
any  data  set.  Table  1  lists  values  of  these  constants  for  monthly  and  annual 
data  for  the  above  regions.  The  coefficient  of  determination,  r  ,  for  the 
various  equations  (implied  in  table  1)  ranged  from  0.50  to  0.99  with  few 
exceptions:  0.2U  fo-  February,  Greenland,  and  0.36  for  March,  Eastern  Siberia. 
(The  value  of  r2  indicates  the  quality  of  fit  between  F  and  N  or  I  and  M: 

1.0  =  excellent  fit,  0  =  no  fit.) 

Stations  were  chosen  in  a  given  area  so  as  to  represent  as  large  a  range 
of  temperature  as  possible.  The  equations  (1)  and  (2)  are  therefore 
applicable  throughout  the  entire  area,  with  only  M  and  N  changing  from  station 
to  station. 
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GREENLAND 


The  constants  _a  and  Jb  may  vary  slightly  with  the  period  of  record.  If 
the  computations  are  carried  out  for  M  and  N  in  degrees  Celsius 
(as  for  table  1),  then  the  value  of  a  equals  160  to  185  days  for  the  annual 
data  or  12  to  18  days  for  the  monthly  data.  In  other  words,  F  (I)  equals  half 
the  days  per  interval  of  time  if  the  minimum  (maximum)  temperatures  is  0°C. 

An  exception  is  for  the  annual  number  of  *rost  days  for  Greenland,  with  a. 
equaling  227  days.  For  the  annual  data,  the  absolute  value  of  b1  or  bo 
appears  to  decrease  with  continentality,  as  from  20  for  Iceland,  a  maritime 
climate,  to  6  for  Eastern  Siberia,  a  highly  continental  climate. 

With  respect  to  the  monthly  data,  the  value  of  a_  tends  to  be  closest  to 
15  for  the  months  of  the  transitional  seasons,  the  spring  or  the  fall,  when 
diurnal  freeze-thaw  cycles  are  most  frequent  in  this  latitude  (see  figure  3); 
usually  the  colder  the  month,  the  higher  the  value  of  _a,  the  warmer  the  month, 
the  lower  the  value  of  a..  If  no  constants  are  given  in  table  1  for  a  given 
month,  then  the  percent  of  frost  or  ice  days  is  either  0  or  100  due  to  the 
relatively  warm  or  cold  temperatures.  In  the  computations  of  the  monthly  data 
the  number  of  frost  (ice  days)  were  limited  to  l-<F(I)^-(n-l),  where  n  =  total 
number  of  davs  per  month. 

A  comparison  between  the  observed  and  the  estimated  annual  number  of 
frost  days  or  of  ice  days  for  test  stations  is  given  in  table  2.  Values  of  r 
for  the  quality  of  fit  between  the  observed  and  estimated  frequencies 
were  >  0.8  for  each  set  of  stations. 


FREEZE-THAW  DAYS 

The  number  of  freeze-thaw  days,  Z,  per  given  interval  of  time 
(information  not  available  in  climatic  summaries)  may  be  found  by  a  variety  of 
methods,  as  by  direct  counting®  or  by  correlation  with  dT. ’*10  in  the  case  of 
Just  the  crossover  of  the  freezing  level,  as  in  this  study,  Z  is  simply  the 
difference  between  the  number  of  frost  days  and  the  number  of  ice  days,  or 


Z  =  F  -  I  (3) 

For  each  area  under  consideration,  the  individual  equations  for  annual  F 
and  I  were  first  obtained  by  means  of  the  appropriate  constants  in  table  1. 
The  annual  Z  may  be  expressed  as  follows: 


r  • 


if 


Alaska 


2  =  1UM  -  ION  -  7 


E.  Siberia 


Z  =  6(M-N)  *-  16 


(5) 


Iceland 


Z  =  18M  -  20N  -  4 


(6) 


Greenland 


Z  =  12M  -  5N  +  59 


(7) 


A  comparison  between  the  observed  and  the  estimated  annual  Z  (table  2) 
shows  somewhat  greater  discrepancies  than  in  the  cases  of  F  and  I,  the  factors 
on  which  Z  depends.  Values  of  r2  which  indicate  the  quality  of  fit  between 
estimated  Z  and  observed  Z  range  from  -^0.60  to  0.76  compared  to  >0.90  for  F 
and  I  (except  for  Iceland  for  which  r2  was  0.84). 

Equation  (5)  implies  that  Z  may  be  correlated  directly  with  (M-N), 
i.e.  AT  for  Eastern  Siberia.  Fraser  (1959)  1  found  a  similar  relationship 
for  Canada,  although  the  cycles  he  investigated  were  for  a  larger  temperature 
span,  -2°  to  1°C,  rather  than  Just  across  the  freezing  level. 


8L.  William.  “Regionalization  of  Freeze-Thaw  Activity."  Annah  of  the  American  Anoclatlon  of  Geographer/,  vol.  14,  pp.  597- 
611,  1964. 

^J.K.  Fraiet.  "Freeze-Thaw  Frequendei  and  Mechanic*!  Weathering  in  Canada,"  Arctic,  vol.  12,  pp.  40-52,  1959. 

10S.S.  Viiher.  "Clime  tic  Mapi  of  Geologic  Interejt."  Bulletin  of  the  Geologic  Society  of  America.  voL  56,  pp.  713-736,  1945. 
^J.K.Fruer.  “Freeze-Thaw  Frequencie*  and  Mechanical  Weathering  in  Canada."  Arctic.  voL  12,  pp.  40-52,  1959. 
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STATION  MODELS  OF  DIURNAL  FREEZE-THAW  CYCLES 


Station  models  of  percent  days  per  month  with  diurnal  freeze  thaw  cycles 
from  January  to  December  are  given  for  selected  stations  in  Greenland,  Alaska, 
Iceland,  the  U.S.S.R,  and  West  Germany  (See  figure  3)*  For  most  of  the  high- 
latitude  stations,  especially  those  close  to  the  Arctic  Circle,  the  peak 
incidence  of  diurnal  freeze-thaw  cycles  occurs  during  the  transitional  seasons 
of  spring  and  fall.  The  stations  are  arranged  so  as  to  show  a  gradual  change 
in  pattern  from  an  extremely  cold  climate  as  Barrow,  Nord,  or  Polar  Station 
where  diurnal  freeze-thaw  cycles  prevail  only  in  summer  to  the  the  relatively 
warm  stations  of  Annette  or  Vladivostok,  which  have  relatively  long  summers 
with  no  freezing  and  the  peak  frequencies  of  freeze-thaw  cycles  in  the 
winter.  The  German  stations  from  Zugspitze  (2962  m)  to  Munich  (532  m)  reflect 
the  effect  of  altitude  on  the  incidence  of  diurnal  freeze-thaw  cycles.  The 
station  model  for  Zugspitze  is  somewhat  similar  to  that  for  Barrow  (or  Polar 
Station).  The  plots  for  the  Icelandic  stations  resemble  that  of  Fichtelberg 
(or  Brocken).  At  Vladivostok,  despite  the  low  latitude,  January  and  December 
are  too  cold  for  freeze-thaw,  whereas,  in  Iceland,  freeze-thaw  occurs  all 
winter. 

Another  set  of  station  models  is  shown  in  figure  4.  This  time  the 
abscissa  is  the  mean  monthly  temperature,  however,  the  ordinate  is  the  same  as 
that  of  figure  3,  namely  the  average  percent  days  per  month  with  freeze- 
thaw.  In  general,  for  a  given  station,  the  daily  freeze-thaw  cycles  per  month 
tend  to  increase  as  the  mean  monthly  temperature  approaches  0°C.  The 
frequencies  per  given  temperature  vary  from  one  station  to  another. 
Nevertheless,  certain  of  the  models  may  sometimes  serve  several  stations  or 
grouos  of  stations. 

Station  models  of  this  type  have  also  been  obtained  for  numerous  other 
stations  in  the  high  and  mid-latitudes.  A  few  of  these  models,  as  well  as  an 
equation  that  provides  curves  that  approximates  some  of  the  data  in  figure  4, 
have  been  given  previously.1^  See  Appendix  for  the  equation.  The  utility  of 
the  models  in  figure  4  is  the  ready  comparison  of  freeze-thaw  regimes  among 
different  stations  and/or  climates. 

The  temperature  limits  for  the  freeze-thaw  regime,  as  well  as  the  peak 
amplitude,  are  climate  dependent.  For  a  highly  continental  station,  as 
Verkhoyansk  (not  shown),  the  freeze-thaw  temperature  regime  extends  from  about 
+  16  to  -16°C  with  an  amolitude  of  about  68  percent  at  about  0°C.  Most 
stations  in  Iceland  have  an  amplitude  of  about  40  to  50  percent  with  a 
positive  temperature  range  of  about  10°C.  Cold  Bay,  Kodiak,  and  Annette  in 
Alaska  have  freeze-thaw  regimes  similar  to  stations  in  Iceland.  On  the  other 
hand,  interior  stations  in  Alaska  have  much  higher  amplitudes,  about  65  to  80 
percent  as  Fairbanks,  McGrath,  and  Gulkana,  with  temperatures  ranges  from  12 
to  -16°C.  For  a  number  of  stations  in  Greenland,  the  plots  (not  given)  show 
much  greater  irregularities  of  pattern  than  for  Iceland  or  Alaska,  possibly 
because  of  the  relatively  short  periods  of  record. 


12R.L.  Vtcxlcr.  "A  Genera]  Climatological  Guide  to  Daily  freezing  Conditions:  Irost  Days,  Ice  Days,  and  Kreezc-Thaw  Days." 
U.S-  Army  Corps  of  Engineers,  t  ort  Belvoir ,  V irginia.  ETL-0287 ,  AD- A 116  771 ,  1982. 
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percent  oays/month  with  freeze-thaw  per  mean  monthly  temp. 


Figure  4.  (continued) 


CONCLUSIONS 


Several  mathematical  and  graphical  models  were  presented  for  estimating 
the  frequencies  of  frost  days,  ice  days,  and  freeze-thaw  days  per  month  or 
year  for  stations  in  Alaska,  Eastern  Siberia,  Iceland,  and  Greenland.  Once 
linear  regression  equations  are  determined  for  the  derivation  of  frost  days 
^and  ioe  days,  respectively,  per  network  of  stations,  these  equations  are  then 
applicable  to  any  site  within  the  area,  given  only  the  mean  daily  minimum  and 
the  mean  daily  maximum  temperatures  for  the  site.  The  frequencies  of  diurnal 
freeze-thaw  cycles  then  may  be  readily  obtained  from  the  difference  between 
the  frequencies  of  the  frost  days  and  the  ice  days. 

Two  types  of  station  models  were  provided,  the  first  the  conventional 
annual  cycle  of  monthly  freeze-thaw  and  the  seoond,  the  same  data  plotted  per 
mean  monthly  temperature.  The  percent  freeze-thaw  for  a  given  temperature 
varies  from  one  station  to  another.  Although  ths  models  tend  to  be 
distinctive,  each  depending  on  lat  tude  (solar  elevation),  altitude, 
continentality,  and  local  conditions,  nevertheless  certain  of  the  models  may 
be  used  to  represent  groups  of  stations. 

The  study  shows  that  diurnal  freeze-thaw  cycles  may  not  necessarily  be 
derived  from  mean  temperatures;  essential  parameters  are  the  mean  daily 
minimum  and  mean  daily  maximum  temperatures.  The  latter  yield  definitive 
information  concerning  frost  days,  ice  days,  and  freeze-thaw  days.  As  a 
consequence,  periglacial  activity  might  be  better  correlated  with  the  mean 
daily  minimum  and/or  mean  maximum  daily  temperatures  per  given  interval  rather 
than  mean  temperature  alone. 
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APPENDIX  A 


A  single  equation  for  the  curves  that  approximates  most  of  the  data  in 
figure  4  is  given  by 

y  =  a  sin  0  (A1 ) 

where  y  =  %  days  per  month  with  freeze-thaw 

a  =  amplitude  of  the  curve  (peak  %  monthly 
freeze-thaw) 

q  =  rb  (T  -  d  +  b)/?b  (Q  in  radians) 

T  =  highest  temperature  associated  with  y  =  0 
H 

Tl  =  lowest  temperature  associated  with  y  =  0 
T  =  temperature  at  amplitude  a 

b  =  b-jor  b2  (see  below) 

bi  =  tH  *  Ta 

b2  =  Ta  -  TL 

d  =  departure  from  0°  C  for  maximum  y 


If  the  curve  is  symmetrical  about  the  peak  percent  frequency  (y  =  a),  then 
b  =  b,  for  the  entire  curve,  otherwise  b  =  b1  and/or  b2  in  turn,  b  is  always 
positive.  For  the  unsymmetrical  plot,  then 

Q  =  ib  (T  -  d  ♦  b^/2^  T-Ti 

6V  =  nr  ( T  -  d  ♦  b2)/2b2  Tf 

Table  A1  lists  the  constants  for  equation  (A1)  for  the  plots  in  figure  4. 
First,  the  plots  were  drawn  manually,  and  the  constants  a,  d,  b1  and  or/b2 
estimated.  The  final  curves  in  figure  4  were  then  generated  by  equation  iAI). 
The  results,  although  somewhat  arbitrary  where  data  points  are  few,  seem  to 
fit  the  coastal  stations  or  the  more  moderate  climates  somewhat  better  than 
interior  stations  as  Fairbanks  or  Gulkana.  The  higher  the  values  of  a.  and  b-, 
or  b2,  the  more  continental  the  climate.  Valdez  appears  to  be  somewhat  of  an 
exception,  with  its  high  a.  of  75  and  low  values  of  b,  and  b2  of  6  and  9.  From 
plots,  as  in  figure  4,  freeze-thaw  regions  may  be  readily  compared  among 
stations  over  the  entire  globe. 


TABLE  A1  Constants  for  Equation  (A1) 


Stations  in  figure  4 


a  d  b'j 


b2 


Barrow 

0.0 

7.5 

Fairbanks 

61 

-2.5 

13.5 

King  Salmon 

55 

0.0 

11.0 

15, 

Grimsstadhir 

45 

2.5 

10.0 

10, 

Natabu 

45 

1.0 

8.0 

Cold  Bay 

50 

0.0 

6.5 

St.  Paul's  Island 

. 50 

..1.0 

_  6.5 

Valdez 

75 

0.0 

6.0 

8 

Anchorage 

60 

-2.0 

10.0 

Summit 

60 

1.0 

9.0 

15, 

Betties 

.  58 

0.0 

_ 12.0 

Gulkana 

70 

-1.0 

12.0 
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